Insulator-Metal transition in the Doped Transition Metal Oxide LaTiO^ 
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The doping induced insulator- metal transition in Lai-^Sr^TiOa is studied using the ab-initio 
LDA-I-DMFT method. Combining the LDA bandstructure for the actual, distorted structure found 
recently with multi-orbital DMFT to treat electronic correlations, we find: (i) ferro-orbital order 
in the Mott insulating state without orbital degeneracy, (ii) a continuous filling induced transition 
to the paramagnetic metal (PM) with x, and (Hi) excellent quantitative agreement with published 
photoemission data for the case of 6% doping. Our results imply that this system can be described 
as a Mott-Hubbard system without orbital (liquid) degeneracy. 



PACS numbers; 71.27.-fa, 71. 15. Mb, 79.60.-i 

The early TMO LaTiO^ is an interesting case of a sys- 
tem exhibiting filling-driven insulator-metal transition. 
Long thought of as a text-book example of a Mott in- 
sulator with antiferromagnetic order, current interest in 
orbital correlations in TMOs has reopened the subject of 
the physics of LaTiO^ [1]. 

According to [2] this material with Ti^~^{3d^) config- 
uration is an electron analog of the cuprates, with G- 
' type AF spin order below Tn — 150 K, and with a re- 
duced spin moment in the AF state. Recent neutron 
scattering measurements reveal an isotropic magnon dis- 
persion, with a small spin gap, = 3.3 meV, putting 
strong constraints on an acceptable mechanism for the 
observed antiferromagnetism [2]. The fact that the sin- 
gle electron in LaTiO^ resides in the t2g orbitals (less 
prone to Jahn- Teller distortions) suggested that the or- 
bital moment might not be fully quenched in this case. 
Attempts to explain the strongly reduced spin moment of 
fi = 0.46/is [2] assuming an unquenched orbital moment 
aligned antiparallel to the spin moment (spin-orbit cou- 
pling in the t2g shell of Ti^^) led to a spin gap much 
larger than that of 3.3 meV observed in INS [2], ex- 
cluding an orbital contribution to the ordered moment. 
Moreover, ferro-orbital order (FOO) has recently been 
suggested [3], in complete accord with the G-type AF 
spin order below T^v, implying lifting of the tig orbital 
degeneracy. 

The metallic phase in Lax-xSr^TiO-i [1] has been 
shown to be a canonical strongly correlated Fermi liq- 
uid [1], putting additional constraints on an acceptable 
model for this system. In particular, this implies a non- 
, degenerate ground state in the paramagnetic metal (PM) 
phase, since any remnant of orbital degeneracy (if it was 
hypothesized to exist) would lead to violation of Landau's 
Fermi liquid hypothesis, and lead to a non-FL metal- 
lic state, in contradiction with observations. AF LRO 
is rapidly suppressed upon Sr doping, and classic sig- 
natures of the Brinkman-Rice (BR) correlated FL are 
seen close to the MIT, suggesting that the AFI state is a 
Mott-Hubbard insulator and that the AFI-PM transition 



is a Mott-Hubbard transition. Given the incomplete un- 
derstanding of the AFI, the mechanism of the MIT still 
constitutes an open problem for theory, in spite of much 
recent activity [4]. 

Two alternatives were proposed to account for the puz- 
zling observations in the AFI phase. Mochizuki et al. [5] 
presented a model for the AF in LaTiOs using a GdFeO^ 
distortion, which, however, was not observed. KhaliuUin 
et al [6] proposed a novel theoretical idea based on an or- 
bital liquid, but the short-ranged orbital RVB-like fluctu- 
ations were not observed either [3]. Moreover, it has been 
shown recently that long range AF order is theoretically 
incompatible with an orbital liquid state [7], deepening 
the controversy regarding the origin of G-type AF order 
in LaTiO^. 

Here, following new work [8] pinpointing the real crys- 
tal structure (RCS) of LaTiO^, we propose a detailed 
theoretical scenario unifying the various experimentally 
observed features in the insulating (I) as well as the corre- 
lated metallic (M) states. Employing the state-of-the-art 
technique LDA-I-DMFT (IPT) (combining the real band- 
structure of LaTiOz with dynamical correlation effects 
of local Coulomb interactions) , we show how a consistent 
description of the doping-induced insulator-metal (I-M) 
transition, as well as very good quantitative agreement 
with photoemission results is obtained. 

The RCS of LaTiO^ (i.e, Pbnm space group) involves a 
structural distortion arising from the tilting of the TiOg 
octahedra around the [110]c axis followed by their ro- 
tation around c, beginning from the cubic perovskite 
structure. The specific distortion of the TiOe octahe- 
dra in LaTiOa is characterized by a large splitting in the 
02 — 02 edge lengths of the octahedral basal plane, giving 
rise to elongation along a and shrinking along b. This im- 
plies FOO and a lifting of tig orbitaldegeneracy. Starting 
from a local picture in the configuration, the occupied 
orbital is constructed as: \a) — 0.77( ^^^J^^^ )±0.636\xy), 
which is non-degenerate [8] . We have carried out a local 
density approximation (LDA) study based on the RCS 
described above. 
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FIG. 1. LDA partial densities of states (DOS) for the real 
crystal structure of LaTiO^. Full line denotes the DOS for the 
non-degenerate ground-state orbital, the dotted line denotes 
the DOS for the higher (almost two-fold degenerate) lying 
t2g orbitals. The dot-dashed line denotes the small eg-DOS 
around Ef- Higher-lying eg and 0-2p bands are not included 
in the Figure. 

To obtain the band DOS, the one-electron scalar rela- 
tivistic Schrodinger equation was solved self-consistently 
within the LDA, using the linear muffin-tin orbitals 
(LMTO) scheme in the atomic sphere approximation, 
with combined correction terms [9]. The exchange- 
correlation potential was parametrized according to the 
von Barth-Hedin scheme [10,11]. Self-consistency was 
reached on performing calculations on a k-mesh of di- 
mension 16 X 16 X 16. The radii of the atoms in the real 
RCS were chosen to minimize their overlap. Values of 
r = 3.7 (La), r = 2.12 (O) and r = 2.75 (Ti) a.u. were 
found to be useful. 

The total LDA density of states, p{uj) = a ~ 
ec((k)), with the LDA (one-electron) Hamiltonian given 
by 

Ho = ^ea(k)cj^^Cka ■ 

Several interesting points are apparent from the LDA 
results. First, the t2g splitting is not negligible, and the 
character of the occupied ground state orbital \a) (not 
shown) is in accord with the result of Ref . [8] . The LDA 
bandwidth of the t2g bands is about 2 eV, and the charge 
transfer energy (LDA) is large, Alda — 7.0 eV, while 
the onsite Hubbard U « 6.0 eV, putting LaTiO^ in the 
Mott-Hubbard category, as expected [1]. Further, the 
Hund coupling, Jh — 1 eV (relevant in the excited (P 
configuration, with the two d electrons on different or- 
bitals, realized upon d — d hopping), and the inter-orbital 
couphng, U' = (U - 2Jh) ^ 4.0 eV. These results (see 



Fig. 1) do not represent an insulating state, this requir- 
ing a proper treatment of correlation effects in the RCS. 
Given the FOO found in LaTiO^, a description in terms 
of a single-band model for the |a) band is untenable, and 
requires a full multi-orbital dynamical mean field calcu- 
lation. To avoid double counting of interactions already 
treated on the average by LDA, wc follow [4] to write, 

Hq = Ho + eia^n^aa , 

with e^„^ = -U{na-a - \) + ^(«q<t - !)■ 
The full many-body Hamiltonian is now given by, 

— Jh ^ ^ia ■ Si/3 ■ 

iafS 

We solve this multiband model in d ~ oo using multi- 
orbital iterated perturbation theory (IPT) [12]. Assum- 
ing no symmetry breaking in the spin/orbital sectors, 

we have Gapaa'{(^) = SapSa-a-'Gacrii^) and Sq/3o-o-'(w) = 

Saf3^crcr'^aa{^)- Further, from the LDA results, the DOS 
for the Bg bands is small near Ep. Taking this together 
with the fact that the O — 2p bands lie 6 eV away from 
Ep^ we keep only the t2g bands in what follows. 

The DMFT solution involves [13] (i) replacing the lat- 
tice model by a self-consistently embedded multi-orbital, 
asymmetric Anderson impurity model, and, (m) the self- 
consistency condition requiring the local impurity Green 
function to be equal to the local GF for the lattice. In 
LaTiOs, in the Sd^ configuration reached by hopping, 
U', Jh scatter electrons between different t2g orbitals, so 
only the total number, nf = „ is conserved in a 

way consistent with Luttinger's theorem. The calculation 
follows the philosophy of the one-orbital IPT, with the 
GFs and the self-energies being matrices in the orbital 
indices. The full set of equations for the multi-orbital 
case are the same as those appearing elsewhere [14], and 
we do not reproduce them here. Local self-energies com- 
puted from this generalized IPT formulation exactly sat- 
isfy the Friedel-Luttinger constraint for arbitrary band- 
fillings, guaranteeing the correct low-energy form of the 
spectral function. The correct form of the spectrum at 
high energy is ensured by the (exact) selfconsistent de- 
termination of the first few moments of the spectral func- 
tion. These equations are solved selfconsistently with the 
LDA DOS as input, until convergence is achieved. 

We now present our results. Figure 2 shows the 
total spectral function for the paramagnetic insulating 
(PI) state of LaTiOs. This does not correspond to the 
true ground state, which is a G-type AFI, but to the PI 
(Mott) state above T/v. The Mott-Hubbard gap is read 
off directly from the DOS as ~ 0.4 eV. Further, the 
renormalized value of the t2g splitting is computed to 
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FIG. 2. LDA+DMFT results for the partial DOS in the 
insulating state, obtained using U — 6.0 eV, Jh ~ 1.0 eV, 
and U' — 4.0 eV in the real bandstructure of LaTiO^. 

be 0.15 eV, implying a reduction in the charge transfer 
energy, Act = i^p — ^d) from 7.0 eV in the LDA to 6.64 eV 
in the correlated solution. This is likely to be an overesti- 
mate, since further downward renormalization of Act will 
result from additional Hartree shifts arising from inclu- 
sion of Madelung {p — d interaction) terms. These results 
constitute a concrete resolution of the questions posed re- 
cently [15], namely, the mechanism for the generation of 
a small-gap Mott-Hubbard insulator in a situation where 
the Coulomb interactions are comparable to the LDA 
bandwidth. In view of our results, dynamical effects of 
coupled spin/orbital correlations in the real structure of 
LaTiOs are of crucial importance for a consistent reso- 
lution of this issue. We also find the occupation of the 
lowest t2g orbital, n^a) = 0.17, slightly higher than that 
of the other two higher- lying (by 0.15 eV) orbitals in the 
t2g sector. 

These results are consistent with having ferro-orbital 
order (FOO) in LaTiOs [3]. Further, since the ground 
state occupied orbital, |a) = ai\xz + yz) + a2\xy) (con- 
sistent with Ref. [8]) is not oriented along one of the 
bonds, the resulting AF exchange is not expected to be 
very anisotropic. A detailed calculation of the Heisen- 
berg superexchange based on the LDA-I-DMFT results, 
along with proper incorporation of spin anisotropies, is 
in progress and will be reported elsewhere. 

To model the paramagnetic metallic (PM) state, 
reached (continuously from the PI in our approach) upon 
doping LaTiOs with Sr, we extended our calculation for 
the PI to: (i) incorporate the change in band-filling along 
the same lines as for generalized IPT for the single band 
case. This implies computing the chemical potential 
shift, as required by Luttinger's theorem, and, (m) con- 
sidering the effects of Sr doping induced static disorder 



(the difference of the on-site potentials on Sr substitution 
is At; ~ 4.0 eV [16]) using the coherent-potential approx- 
imation (CPA) combined selfconsistently [17] with the 
multi-orbital IPT described above. Furthermore, to ob- 
tain quantitative consistency with spectroscopic results, 
we found it necessary to include the (small near Ep) DOS 
of the Cg bands foimd in the LDA (Fig. 1), as well (see 
below) . 

Generically, hole doping the (Mott) PI creates a very 
sharp, quasicoherent Fermi liquid peak at the renormal- 
ized Ep, with dynamical spectral weight transfer over 
a large energy scale from high to low energy (caused by 
U, U', Jh within DMFT), driving a first-order Mott tran- 
sition in the generic case. In this situation, static, poten- 
tial disorder is expected to lead to a moderate broadening 
of the sharp FL peak, with spectral weight transfer from 
low to high energy [17], and, more importantly, to con- 
vert the first-order Mott transition into a continuous one. 
Thus, the character of the MIT, as well as details of the 
actual one-particle spectrum depend sensitively on how 
the competition between Mott-Hubbard localization and 
itinerance in a multi-band system is affected by doping 
induced static disorder. 

With this in mind, we show the photoemission (PES) 
and inverse-PES (IPES) spectra for the lightly (6 = 
x+y = 0.06) doped, PM state of La^Sri^^TiO:i+yi2 [18] , 
along with the experimental result for Ippsi^) in Fig- 3. 
Clearly, excellent quantitative agreement between theory 
and experiment is found over the entire range —2<uj< 
Ep. In particular, the lower Hubbard band, the quasi- 
coherent FL feature, as well as the details of the lincshape 
are all in excellent agreement with experiment. In the 
IPES spectrum, we are able to resolve the a, (3 peaks 



0.4 




FIG. 3. The photoemission (PES) and inverse photoemis- 
sion (IPES; in inset) spectra for both, pure (dashed) and 
doped (solid) LaTiOz ■ Excellent quantitative agreement with 
experimental results (PES: squares) of Ref. [18] is clear. 
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alongwith the much greater signal intensity and the cor- 
rect detailed form of the lineshape [19]. Thus, our results 
represent a quantitatively accurate representation of the 
complete one-electron spectral function in the PM phase 
of Lai-xSrxTiO^j^y/2- Lastly, from the Fermi liquid self- 
energies (not shown), the effective mass renormalization 
is estimated to be m*/mt, = 1.5 for 5 = 0.06 in good 
quantitative agreement with the frequency-dependent 
renormalization estimated from PES [18]. Actually, the 
total renormalization includes additional renormalization 
(factor of 2) coming from the k-dependence of the self- 
energies. On general grounds, we do expect such renor- 
malization in a case of proximity to the AF transition; 
theoretically, an extension to include non-local correla- 
tions is required to address this issue. We have not done 
this here. 

We stress essential differences between our work and 
earlier works. In Ref. [4,5], an incorrect crystal structure 
was used for the LDA calculation. Moreover, the strong 
Jh — 1.0 eV along with the inter-orbital U' seem to have 
been completely ignored in their correlation (DMFT) cal- 
culations. Turning to Ref. [6], we believe, given the 
fact that a multitude of experiments reveal correlated 
(Brinkman-Rice) FL responses (PM), that the basic Lan- 
dau hypotheses force one to abandon a picture based 
on orbital degeneracy in the PM phase as well. We re- 
emphasize that AF-LRO has recently been shown to be 
impossible [7] for a id} system with triple orbital {t2g) 
degeneracy. Very recently, Pavarini et al. [21] have stud- 
ied a host of Zd} systems within LDA+DMFT, obtain- 
ing similar conclusions for the Mott insulating state in 
LaTiOz- However, the S'r-doping induced I-M transition 
has not been studied there. Additionally, using QMC to 
solve the DMFT equations forces a restriction to high- 
T = 770 K [21]. Given the strong T dependence of the 
correlated Fermi liquid scale within DMFT, as well as the 
much lower FL coherence scale observed experimentally 
in doped LaTiO^ [1,13], (in contrast to QMC results), we 
believe that a quantitative comparison with experiment 
at low T requires the solution of the DMFT equations at 
lower T ~ 150 A', as we have done. Finally, we draw at- 
tention to the excellent quantitative agreement between 
theory and experiment for the one-electron spectral func- 
tion, obtained using LDA-|- (multi-orbital) DMFT. To 
our knowledge, this is the first calculation achieving such 
unprecedented agreement between experiment and the- 
ory for this system. 

To conclude, we have studied the doping induced Mott 
transition in the system Lai-xSvxTiOz using the 
ab-initio LDA-I-DMFT method. Starting from the real 
crystal structure (with lifting of t2g orbital degeneracy), 
we have explicitly demonstrated the nature of the carrier 
driven MIT in this system. Our results support ferro- 
orbital order in the insulating state, consistent with re- 
cent suggestions [3], and with the observed G-type AF 
order. In conjunction with [8, 21], these results favor a 



description of the physics of LaTiOy, in a picture of non- 
degenerate t2g orbitals. Within this picture, in the PM 
state, excellent quantitative agreement with published 
PES and IPES results is obtained, providing a complete 
description of the evolution of one-particle excitations 
across the MIT. Similar ideas can be fruitfully applied to 
study the doping-driven MI transitions in other TMO- 
based systems [20] with coupled spin-orbital correlations. 
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